Experimental characterization of vibrated granular rings 
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We report an experimental study of the statistical properties of vibrated granular rings. In this 
system, a linked rod and bead metallic chain in the form of a ring is coUisionally excited by a 
vertically oscillating plate. The dynamics are driven primarily by inelastic bead-plate collisions and 
are simultaneously constrained by the rings' physical connectedness. By imaging many instances 
of the ring configurations, we measure the ensemble averages and distributions of several physical 
characteristics on the scale of individual beads and composite ring. We study local properties such 
as inter-bead seperation and inter-bonds angles, and global properties such as the radius of gyration 
and center of mass motion. We characterize scaling with respect to the size of the chain. 

PACS numbers: 05.40. a, 81.05.Rm, 82.35. Lr 
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I. INTRODUCTION 

In the last decade the study of granular materials has 
become a mainstream pursuit amongst physicists who 
have summarized the advances in recent reviews Q, |^ |^ . 
The broad spectrum of static and dynamic behaviors of 
granular media are being studied by a combination of 
state-of-the-art experiments, computer simulations and 
theoretical approaches. Shaken granular materials are 
one example of a comprehensively studied subset of this 
spectrum. Experiments and simulations on shaken grains 
have demonstrated a variety of phenomena such as local- 
ized structures , extensive patterns 0| , clustering 
B. Hi , c ompaction 0, 0| and non-Maxwellian statistics 

Granular chains consist of hollow spheres connected 
by rods, i.e., grains constrained by a chain backbone. 
In granular chains we get the juxtaposition of the hard 
grains of granular physics with the flexibility of polymers. 
The chains can be excited, say, by collisions with a vi- 
brating plate and can be observed visually. As such they 
afford a simple system for the direct measurements of 
micro- and macro-scopic variables. Experiments on gran- 
ular chains, supported by theoretical models and simu- 
lations, have probed the diffusive relaxation of a simply 
knotted chain |l3j . the dynamical behavior of a verti- 
cally excited hanging chain the entropic tightening 
of a once-twisted "fi gure -S" ring and the spontaneous 
formation of spirals [iq . 

In this paper we report our study of vibrated granular 
rings. Our experimental system consists of a vertically 
oscillated plate on which we place a granular ring con- 
structed from a ball chain; a schematic is shown in Fig.^ 
Rings are constructed by joining the two free ends of a 
chain. A chain with free ends suffers unequal excitation 
along its length since beads near the ends are less con- 
strained and thus preferentially excited. Rings, which 
are devoid of free ends, restore the equality in excitation 
amongst the beads, and it is for this reason that they are 
particularly attractive to study. 

In our experiments we control the amplitude and fre- 
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FIG. 1: A schematic of our experimental apparatus. 



quency of the harmonically oscillated plate. A chain ini- 
tially resting on the plate becomes airborne when the 
plate imparts to it an acceleration that exceeds the grav- 
itational acceleration. Thereafter, collisions between the 
plate and chain appear to randomize the motion and 
shape of the granular ring. The collisions of different 
parts of the chain are asynchronous with the phase of the 
plate's oscillation. Further collisions amongst the beads 
on the chain occur frequently. 

We observe the ring from above at a random phase 
of the oscillation. We acquire images which are two- 
dimensional projections of the position of the ring. From 
these images we obtain the spatial coordinates and order 
of the beads in the ring. Since the images are sequential 
and time stamped, we are also able to calculate differ- 
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ences and therefore dynamics. We make a large number 
of position measurements for rings under vertical vibra- 
tion. From this data we study the geometrical proper- 
ties on the bead scale, i.e., local characteristics, and on 
the scale of the chain, i.e., global characteristics. The 
local properties include the inter-bead separations, the 
bending angle between two consecutive rods, the aver- 
age curvature and the dynamics of the beads. For the 
global characteristics we focus on descriptors of the size 
and overall motion of the ring. These are the radius of 
gyration and the dynamics of the rings' center of mass. 

This paper is organized as follows. In Section ITU we 
describe the experimental apparatus, the specific design 
considerations, the image acquisition and analysis proce- 
dure. Our experimental results depend at most on the 
ordered positions of the beads in a ring and on the data 
acquisition rate. We have divided the results into local 
and global properties which are discussed in Sections IIIII 
and IIVI respectively. Local (global) properties are rel- 
evant to the scale of the bead (ring). We discuss the 
the inter-bead separation, the inter-rod angle, the aver- 
age curvature and the bead dynamics as local properties. 
For global characteristics we discuss the radius of gyra- 
tion and the average motion of the ring center of mass. 
A conclusion and an outline of future directions follows 
m Section El 



II. EXPERIMENT 

In this Section we describe our experimental apparatus 
and our data extraction procedure. In overall design, the 
apparatus is comprised of an electromechanical shaker 
that drives a plate which is observed by a CCD camera. 
A schematic of our experimental setup is shown in Fig. ^ 



A. The electromechanical vibration assembly 

A commercial electromechanical shaker driven by a si- 
nusoidal voltage input is used to oscillate a plate. The 
shaker operates electromagnetically with the current- 
carrying coil armature oscillating in the strong magnetic 
field of permanent ceramic magnets. The electromechan- 
ical shaker is coupled to the plate by a shaft but is other- 
wise vibrationally isolated from the rest of the apparatus. 
A massive aluminum platform that is rigidly bolted to 
the laboratory floor acts as a low-pass filter of vibrational 
noise. The shaft is guided by the bore of a precise square- 
section linear air bearing. In addition to linear motion, 
the square-section shaft decouples the plate from the ro- 
tational torque in the drive. The driving is thus accu- 
rately uni-axial and a leveling plate facilitates alignment 
with gravity. The shaft is coupled to the plate by two 
concentric cylindrical tubes. A sensitive accelerometer is 
attached centrally on the underside of the plate. 

The plate is driven harmonically with amplitude A and 
angular velocity oj — 27r/, where / is the driving fre- 



quency. The dimensionless peak acceleration V = AuP' /g, 
where g is the gravitational acceleration, and the fre- 
quency / specify the driving conditions for the experi- 
ment. The shaker is operated at / = 16 Hz and T = 2.12 
for all the experiments described in this paper. The 
shaker is protected from bearing the combined weight 
of the plate and shaft by adjustable counter-weighting 
springs such that the drive oscillates about its natural 
equilibrium position. 



B. The plate and the granular rings 

The plate is circular with a diameter of 11.50 inches 
and is made of stress-relieved aluminum. The primary 
design requirements are light weight, flexural stiffness 
and surface hardness. Aluminum's low density and com- 
parable flexural properties to most metals make it a suit- 
able choice. On oscillation we expect that bending modes 
of the plate are excited. The higher the stiffness the 
higher the frequecies excited. Since our experiments have 
drive frequency / ~ 10 Hz, we designed the plate so that 
its excitable bending modes had much higher frequencies. 
This is accomplished by tapering the lower surface of the 
plate and by boring out cylinders from this surface in a 
space-filling azimuthally symmetric manner. The cylin- 
ders bored out of the aluminum have various depths and 
diameters. The plate is at its thickest 1.50 inches decreas- 
ing to 0.25 inches at the boundary. These considerations 
lead to a selective damping of pure modes of the shaken 
plate. A further safeguard against exciting modes of the 
plate is achieved by forcing the plate at two annular re- 
gions rather than at the center or at the edge. The an- 
nular regions are at radial locations 1.25 < r < 2.00 and 
4.88 < r < 6.00 inches. Surface hardness needed to with- 
stand the many metal on metal collisions of the chains 
hitting the plate is attained by hard black anodization of 
the aluminum plate. 

The granular chains are commercially available ball 
chains. They consist of hollow, approximately-spherical 
shells that are connected by dumb-bell shaped rods. We 
use two types of ball chains in our experiments. The 
smaller is made from nickel-plated stainless steel rods 
and beads. The beads have a nominal diameter of ^ 
inches. The larger is a ^ inch diameter brass chain. Thus 
there are typically 100 beads across the plate diameter. 
The chains are connected end to end to make rings. The 
splicing together of the free ends is natural in the sense 
that a rod from one of the end beads is inserted into the 
shell of the other end bead. The joined free ends could 
seldom be identified in the constructed ring. The tightest 
rings that we could make have 8 beads. Hence the an- 
gle between two consecutive rods is constrained so that 
< 6* < 50°. Our experiments use rings with N beads 
where 13 < < 309. The rings are laterally confined to 
the plate by an acryllic boundary. 
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C. Data extraction 




(a) 





Our data consist of digital images of the granular ring 
on the plate. We illuminate the chain using approxi- 
mately normally-incident light from far above the plate. 
To minimize reflection from the background, we paint 
the acryllic boundary black. Additionally, before black- 
anodizing the flat surface of the plate is roughened by 
sand blasting. The very shallow surface inhomogeni- 
ety helps randomize even simple collisions with the plate 
and diffusively scatters far- field light. The poorly reflec- 
tive background reduces noise in the digital image. The 
normal far-field illumination gives an almost uniform re- 
flected light intensity across the plate area. Furthermore, 
the light reflected by the beads is primarily from the 
center of the beads. Non-normal lighting preferentially 
illuminates crescent-shaped areas on the surface of the 
beads. A representative raw image is shown in Fig.|2jL. 

Typically, our images are acquired at sampling rates 
of 0.25 — 16 Hz and have a resolution of approximately 
1000 X 1000 pixels. Depending on the data set we have 
between 8.8 and 13.0 pixels to the bead diameter d. The 
reflected light is brightest from the the center of the 
beads and diminishes radially outwards. A sequence of 
standard image analysis procedures such as background 
subtraction, intensity rescaling, high-pass thresholding 
and smoothing result in reducing the data into roughly 
N X 7 X 7 non-zero pixels for N beads. Local maxima 
in the intensity field are then determined and ranked by 
their peak intensity. The intensity field of each set of 
5x5 pixels centered about a local maximum is fit to an 
azimuthally symmetric Gaussian. The fit parameters of 
the position of the Gaussian peak are taken as the sub- 
pixel coordinates of the bead. The fit to the standard 
deviation of the Gaussian gives an estimate to the posi- 
tional accuracy of about 0.05 bead diameters. Figure |2; 
shows a blow up of part of the raw image in Fig. and 
our determination of the bead centers. 

In most of the images, our data extraction procedure 
resulted in more than N local maxima. Most of these 
occur when a rod appears brighter than a bead. A rod 
between two beads is visible in Fig. [2;. The frequency of 
this occurrence is about 1% of the number of beads. Since 
the distance between the position of a rod and an adja- 
cent bead is smaller than the minimum distance between 
two adjacent balls, we eliminated these by iteratively ap- 
plying a strict minimum distance criterion. 

The next step is to order the N coordinates so that the 
ring can be reconstructed. The ordering algorithm uses 
the geometrical restrictions imposed by the interconnect- 
ing rods. Given two adjacent beads, the third is searched 
for in the sector defined by the maximum inter-bead sep- 
aration and angle. This efficient greedy algorithm re- 
quires N'^ operations to order an N bead chain. The 
overall accuracy of our data extraction process exceeds 
96% ,i.e., in less than 4% of the images are we unable to 
reconstruct the chain. There are various reasons for this. 
The beads are constructed by crimping sheet metal and 



FIG. 2: (a) Raw data image, (b) reconstruction by identifying 
the bead coordinates (in pixels) and order, and (c) overlay of 
a section of the raw image and the position data (•) for the 
center of the beads. Locally bright regions between beads in 
(c) correspond to reflection from an interconnecting rod. 
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so have a seam running along one half of a great circle. 
The seam is comparatively duller in reflected light than 
any other part of the beads and so in a few images our 
extraction process misses these beads. Since the images 
are two-dimensional projections of the chain, the strict 
inter-bead distance criteria can appear to be violated. In 
these cases our algorithms may allow a false bead posi- 
tion. In Fig. I^h, we show a complete reconstruction of 
a chain from the raw image data in Fig. EK- Given the 
small number of misidentifications, we are confident that 
these issues have a small effect on the results reported 
below. 



III. LOCAL PROPERTIES 

In this Section we present our results on the local prop- 
erties of vertically vibrated granular rings. Local prop- 
erties are relevant to the scale of the bead and here we 
discuss the inter-bead separation, the inter-rod angle, the 
average curvature and the bead dynamics. We noted ear- 
lier that the chains are constructed out of hollow shells 
that are connected by rods. This connection imposes 
constraints on the minimum and maximum separation 
between beads, on the maximum bending that 3 adja- 
cent beads can sustain, and on the displacements that 
occur on excitation. 



A. Inter-bead separations 

The inter-bead separation 6 = |ri+i — r^j is measured 
between every two consecutive bead positions for brass 
and nickel rings, ri is position of bead i in two dimen- 
sions. Hence, there are N measurements of 5 for an A^- 
bead ring. Our data is based on no fewer than 5200 mea- 
surements of 6 ior N = 13 and more than 10^ samples for 
N > 256. From these data, we construct the probability 
distribution function of the inter-bead separations PiS). 

In Fig. 131 we plot the mean (S) and standard deviation 
as as a function of N. For N sufficiently large, {S) and ag 
are independent of N. The distributions P{S) are similar 
in gross feature for the brass and nickel rings for the var- 
ious lengths N. Combining the data for 30 < iV < 309, 
we obtain the distribution function shown in Fig. ^ The 
distribution has a sharp cut-off at 5 1.4 corresponding 
to the maximal extension between the beads. A much 
softer small 5 cut-off appears at 5 w 0.8. The minimal 
distance in an accurately two-dimensional chain would be 
at S = 1. The quasi- two-dimensional motion of the chain 
and our measuring the motion on a planar projection al- 
lows values of (5 < 1. The peak of the distribution at 
small separations is possibly due to bead-bead frictional 
effects. Between the small-5 peak and the upper-(5 cut-off 
we find that the distribution is roughly flat. These inter- 
mediate separations occur with signiflcant and approxi- 
mately equal probability, perhaps because of the constant 
rod-bead friction at these extensions. 
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FIG. 4: The probability distribution P{5) for the inter-bead 
separations. The distribution combines data for rings with 
lengths 30 < Af < 309. The dashed line S = 1(1.4) in the 
theoretical lower (upper) bound. 



B. Bending angles 

Given the positional ordering r^-i , r^, r^+i of three con- 
secutive beads in the ring, we describe the bending by 



cos 



\dr,^i\\dri 



(1) 



where dr^ = ri+i — r^. There are N values of cos 9 for 
every A^-bead ring. As with the inter-bead separations, 
we have approximately 5200 measurements of cos 9 for 
iV = 13 and more than 10^ samples for N > 256 from 
which we determine the distribution function. In Fig. [Sj 
we plot the mean (cos 9) and standard deviation acose 
of the bending distribution functions for granular rings 
as a function of N. As N increases, (cos 6*) — > 1 and 



Ccose 0. This large N behavior imphes that a long ring 
is primarily composed of locally straight segments. For N 
approaching the smallest ring size of 8 beads, (cos 9) 
0.7 consistent with the tightest bending. 

The large N limiting behavior of the mean bending an- 
gle is particularly interesting. In Fig. we plot the dif- 
ference between the limiting value of unity and the mea- 
sured (cos 9). For an ideal circular ring with N beads, 
1 — (cos 9) ^ N~'^. We find that although the granular 
ring becomes locally straight as its length increases, it 
does so very differently from a circular ring. Our data 
show a gentler approach with 



1 - (cos 61) - N' 



-0.6 



(2) 



This implies that the chain curvature decreases slowly as 
the chain length increases. 

The bending distributions for representative small and 
large rings are shown in Fig. [7| The distributions peak 
near cos 0^1 and monotonically decrease for smaller 
cos 9. The peak and range depends weakly on the num- 
ber of beads N. The peak probability for bending angles 
is higher for the longer rings whereas the range of bend- 
ing angles is larger for the shorter chains. There are 
significant quantitative differences (not shown) between 
the distributions for brass and nickel rings but the qual- 
itative features are similar. The principal difference is 
at large angles where the greater flexibility of the brass 
chains increases the range of variation in cos 9. 



C. Average curvature 

The vibrated rings sample many microscopic configu- 
rations subject only to the stretching and bending con- 
straints imposed by the interconnecting rods and by mu- 
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tual self-exclusion. One of the properties of the configu- 
ration is the overall extent to which the ring is curved. To 
probe this feature we use the ordered positions to com- 
pute the average squared curvature K which we define 
as 



K = 



'k{s)-k{s) ds 



(3) 



kcirclel-^) ' kcircle('5) ds 

where k is given by 

k(s) = ^; ris)^ix{s),y{s)). (4) 

In Eq. 01 s is the arc length along the chain and r(s) 
is a unit-speed reparametrization of the chain. Then, 
by definition \dr/ds\ = 1 so k(s) measures the way 
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FIG. 7: The probability distribution of the bending angle 
P{cos 6) for rings with lengths iV = 32(o) and N = 128(n). 
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FIG. 8: The mean {K) (o) and standard deviation gk (•) of 
the squared curvature, K vs ring length N . The horizontal 
dashed line K = 1 represents a circle. The solid lines are 
proportional to A''^'^ and A''^. 

the chain turns, kcircie(s) is obtained by computing 
Eq. ^ for a circle with circumference equal to the arc 
length of the chain. The local squared curvature is 
k{s) = k(s) • k(s). Averaging over the ring we get 
§ k{s) ds/ §ds. Similarly for the circle with circum- 
ference equal to the chain's length, we have an average 
squared curvature § kdrdeis) ds/ § ds. Normalizing the 
average squared curvature of the ring by that of an equal 
arc length circle, we get the quantity K defined in Eq.|31 
As a procedure, for each ring, we first calculate its cir- 
cumference by adding the separations between adjacent 
beads. We then reparameterize the positions of the beads 
by the arc length s so that = {xi,yi) = {x{si),y{si)) 
such that the ring has unit circumference and is traversed 
with unit speed. Since /kcirc/e(s) • ^circie{s) ds = Att'^ 
for a circle of unit length and speed, we have normal- 
ized K such that as the ring approaches a circle K ^ 1. 
Hence K measures the departure from a circular shape. 
In practice, the integral in Eq. |31 is replaced by a sum- 
mation over the N discrete bead positions and k(si) is 
approximated by the difference 

, , ^ rf^ri dvi+i - dvi 

= ^ = TT- 7\-TT: ^ -y (5) 

Using in excess of 6800 images or approximately 425 
realizations per chain, we have measured the mean and 
rms values of K for brass and nickel rings for several N . 
In Fig. |S1 we plot the mean {K) and standard deviation 
ax as a function of N . The two smallest rings have 
= 13 and 20. They are approximately circular and 
tightly constrained, resulting in (A') « 1; the smallest 
ring has N ~ 8 and is practically rigid. For larger A^ we 
find that the scaling exponent /i ~ 2 in the power-law 
(AT) ~ A^ for 20 < A^ < 256. This imphes, in agreement 
with Eq. 2, that as A^ increases, the ring increasingly 



deviates from a circle. The standard deviation aK from 
the data shows a scaling relation ax ~ A^^ with A ~ 1.5. 
This observation is consistent with Eq. |21 

D. Bead displacements 

The dynamics of the beads are studied by measuring 
the displacement of a bead as a function of time. It is well 
known that the mean square displacement of a randomly 
excited particle scales linearly with time while ballistic 
motion results in quadratic scaling. It is interesting to 
probe how a bead constrained to the ring behaves. 

We restrict our analysis to images acquired at 16 Hz 
which is sufficiently fast so that bead positions between 
two consecutive images show small changes. Conse- 
quently, we can identify beads in an image with the pre- 
ceding image. Discerning the bead identity is increas- 
ingly more difficult with slower image acquisition rates 
since the bead displacements are quite large. The or- 
dered position data of the beads in the ring is obtained 
as described earlier in Section FlI CI Then assuming the 
bead identity in the first image, we order the bead posi- 
tions of the second image so that they correspond to the 
preceding image. Since there is only 1 plate oscillation 
and therefore on average 1 chain excitation between con- 
secutive images, we are confident that we have accurately 
identified the beads. Additionally, we visually verify the 
ordering for a subset of the data sets. We implement 
our bead identity algorithm by translating the data from 
both images to a common center of mass and then min- 
imizing the mean square position difference between the 
two images up to cyclic permutations and orientation of 
the second image. The procedure is iterated until all the 
images are properly ordered. 

The squared displacement that bead z in a ring 
with A^ beads undergoes between time t and t + t is 
Arf{t,T,N) = (ri(r + t)-ri(T))-(ri(r + t)-ri(T)). For 
the data, the times r and t are chosen to coincide with 
the image acquisitions, i.e., multiples of 1/16 s, and the 
displacement is measured in units of the bead diameter. 
Note that both images are translated to the origin, i.e., 
the center of mass is subtracted from the bead coordi- 
nates before calculating the squared displacement. 

Averaging over all beads i in the ring and over possible 
values of t, we define the mean squared bead displace- 
ment: 

{Ar^){t,N)^{ArUt,T,N))r,,. (6) 

We probe the time dependence of (Ar^) by assuming 
a length dependent constant factor l3{N) such that the 
mean squared displacements are the same for the short- 
est experimentally observed time t = 1/16 s. In Fig. |51 
we plot {Ar^) / (3{N) versus t. For times corresponding to 
about 10 plate oscillations, (Ar^)//?(A) increases linearly 
with time for all rings. This behaviour is consistent with 
diffusive motion. Thus, although the ring is mechan- 
ically and athermally driven, collective properties may 
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FIG. 9: The mean squared displacements (Ar^) of the beads 
as a function of time t for various ring length A^. The displace- 
ments are scaled by a ring length-dependent constant P{N) 
so that all (Ar^) are the same for the shortest experimentally 
observed time t — 1/16 s. The dashed line is proportional to 
t. Th 
solid 
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FIG. 10: The probability distribution function for the 
Ar^/(Ar^> at t = 1/16 s for all ring lengths N. The dashed 
line has slope —4. 

still behave as if they are randomly (thermally) forced 
[T^ ITsj l . For shorter times the mean squared displace- 
ments increase faster than linearly with time for all but 
the smallest ring which is consistent with the expected 
crossover from very short time ballistic motion to long 
time diffusive dynamics. 

For the shortest observation interval of t — 1/16 s, 
we study the distribution of the squared displacements. 
Since each bead in the ring collides with the plate, on 
average, once per oscillation, it is highly probable that 
several beads in the chain will have altered direction due 
to a collision. Thus, it is unlikely that the 1/16 s interval 



is sufficiently short to deduce an instantaneous velocity. 
Since the bead-plate collisions occur over a range of phase 
in the oscillation, however, the beads encounter a range 
of displacements. Large displacements correspond to col- 
lisions and smaller displacements to beads being dragged 
along by the chain connectedness. In Fig. ^| we plot 
the probability distribution function of Ar^ measured at 
t — 1/16 s. The squared displacement is normalized 
by its mean value. Our data show that small squared 
displacements are highly probable decaying slowly until 
0((Ar^)). For larger displacements the probability de- 
cays more rapidly and the large displacement tail decay 
is consistent with (Ar^)~^. 

IV. GLOBAL PROPERTIES 

Global characteristics of granular rings probe the over- 
all conformation and dynamics of the chain. Below we de- 
scribe the conformation, roughly speaking, by the chain's 
size. The size, defined to be the mean radius of gyration, 
is obtained by averaging over many observations. The 
global dynamics are described by the motion of the cen- 
ter of mass of the ring. 

A. Radius of Gyration 

A measure of the global size of a long chain or a poly- 
mer is its radius of gyration. Since the chain's configura- 
tion is constantly changing owing to collisional excitation 
with the plate, better measures are the statistical aver- 
ages of radius of gyration. Given the N position coordi- 
nates Ti = {xi, Hi) for the beads i — 1, .., N, we calculate 
the radius of gyration Rg from 



i=l 

where rem = (xcm^ycm) is given by 



rcm^^^r;. (8) 

1=1 

Since our coordinates are in units of the bead diam- 
eter d, the radius of gyration Rg is in these units. We 
have measured Rg for granular rings for several values of 
N. Since there are no self-intersections, the ring can be 
likened to a self avoiding walk that returns to itself. 

Averaging over approximately 450 realizations per 
granular ring of length iV, we measure the mean and 
rms values of Rg. In Fig. ^2 we plot {Rg) and as a 
function of TV for 13 < iV < 309. In these experiments 
the lateral boundary has a radius of 53.2 bead diameters. 
This is an absolute upper bound on (Rg) and is shown by 
the horizontal line in Fig. ^2 Our (Rg) data for nickel 



FIG. 11: The mean (o) and standard deviation (•) of the 
radius of gyration vs the length of the rings N. The dashed 
horizontal line is the radius of the plate. The two solid lines 
have slopes of 3/4 (top) and 1 (bottom), respectively. 
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FIG. 12: The mean (o) and standard deviation (•) of the 
center of mass velocity. The straight lines are power-law fits 
to the data with slopes -0.67 (mean) and -0.59 (standard de- 
viation), respectively. 



and brass chains are consistent and increase with N. A 
power-law fit for the range 20 < iV < 219 gives an ex- 
ponent ly = 0.85 ± 0.05. The a^^ data in Fig. II II shows 
a linear increase with iV for 20 < < 219 with a sharp 
drop for larger N. The sharp drop is associated with the 
constraining effect of the lateral boundary. 

Self avoiding walks have universal scaling characteris- 
tics in the large-A^ limit. In particular, ( Rf, ) ~ with 
Flory exponent ly — j in two dimensions [13, ■ Exper- 
iments with surface adsorbed DNAT^I, polymer molecu- 
lar monolayers 19] and connected plastic spheres 20] have 
measured exponents consistent with the 2D Flory value 
u = J. Both the DNA and plastic chain had free ends 
and by being constrained to 2D were self-avoiding. Our 
measured exponent v = 0.85 ± 0.05 is for a ring (likened 
to self avoiding walk that closes on itself) and is some- 
what higher though not inconsistent with these experi- 
ments. The fjE. increase linearly with N in the range 
20 < iV < 219. 



B. Center of Mass Motion 

Under oscillation the rings sample many configurations 
with different positions of their center of mass. Since 
the plate is accurately level, the rings do not drift grav- 
itationally. Any motion of the center of mass is then 
associated with the overall mobility of the rings. We 
study the motion of the center of mass of nickel rings 
with 13 < < 309. The data are acquired at 16 Hz 
which is sufficiently fast for probing the motion of the 
ring's center of mass. Given the center of mass we 
calculate the magnitude of the average velocity v given 



by 

" = V^'^- 

^From approximately 400 images per ring we measure 
the mean and rms values of v. In Fig. 1121 we plot the 
mean (u) and standard deviation for several N . For 
our largest ring N = 309 significant interaction with the 
lateral boundary appears to be responsible for the un- 
expectedly large {v) and a^. Assuming scaling relations 
{v) - TV and ct„ ~ N'^ power-law fits for 13 < iV < 256 
give ?7 = -0.67 ± 0.08 and C = -0.59 ± 0.14. The speed 
V varies from just over 1 bead diameter/s for = 13 
to about 1/10 bead diameters/s for N = 256. As N 
increases more and more segments of the ring move inde- 
pendently resulting in large cancelations to contributions 
to center of mass motion. Qualitatively, this diminished 
mobility with increasing chain length resembles the be- 
havior observed for a random polymer where the diffu- 
sivity is inversely proportional to the molecular weight 

m 

As with the bead displacements studied in Sec- 
tion llll D1 we can study the dynamics of the ring center of 
mass. The squared displacement that ring center-of-mass 
undergoes between time r and r -I- i is A_R^^(i, r, N) — 
(rcm(T + t) - rcm(T)) • (rcm(T + t) - rcm(T)). For the 
data, the times r and t are chosen to coincide with the 
image acquisitions, i.e., multiples of 1/16 s and the dis- 
placement is measured in units of the bead diameter. Av- 
eraging over all possible values of t we define the mean 
squared ring center of mass displacement: 

{l^Rl^){t,N) = {^Rl^{t,T,N))r. (10) 

The time dependence of (AR'^^^) is plotted in Fig. ^| 
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FIG. 13: The mean squared displacements (Ai?^m) of the 
ring center of mass as a function of time t for various ring 
length A'^. The displacements are scaled by an ring length- 
dependent constant a{N) so that all (Aiicm) are the same 
for the shortest experimentally observed time t = 1/16 s. The 
dashed lines are proportional to t and t^. The (AJJ^m) from 
the three longest rings N — 219, 256 and 309 are shown in 
magenta, cyan and maroon open square symbols respectively. 



Here, we have assumed a length dependent constant fac- 
tor a{N) such that the mean squared displacements are 
the same for the shortest experimentally observed time 
t = 1/16 s. 

For times corresponding to less than about 10 plate os- 
cillations, the {ARl^)/a{N) increase quadratically with 
time for all rings except the three largest rings with N — 
219, 256 and 309. The constraining effect of the plate's 
lateral boundary explains the less than quadratic time 
dependence of {/S.R^^) / a{N) for large rings. For times 
in excess of about 10 plate oscillations, the {AR^^) / a{N) 
approaches a linear increase with time. The short time 
behaviour is consistent with ballistic motion of the ring 
center of mass. This suggests that at short times the 
motion of the ring center of mass is dominated by coher- 
ent large bead displacements. The longer time behaviour 
averages over several, apparently random, large displace- 
ments resulting in the approach to diffusive motion. 

The distribution of the squared displacements of the 
ring center of mass is studied for the shortest observa- 
tion interval oft = l/16s. At this time scale the cen- 
ter of mass motion is ballistic and so an average veloc- 
ity proportional to the displacement may be infered. In 
Fig. 1141 we plot the probability distribution function of 
Ai?^„ measured at t = 1/16 s. The squared displacement 
is normalized by its mean value. Our data show that 
small squared displacements are highly probable with the 
probability decaying consistent with (Ai?^„)^^-^. The 
small squared displacements of the ring center of mass 
arise from averaging the almost independent motions of 
the beads on the ring. The less probable large squared 
displacements occur from rarer coherent motions of seg- 



FIG. 14: The probability distribution function for the 
AR1:m/{ARIm) at t = 1/16 s for all ring lengths N. The 
dashed line has slope —1.5. 



ments of beads on the rings. 



V. CONCLUSION 

In this paper we presented an experimental study of 
vibrated granular rings. Our study elucidates the statis- 
tical characterization of various local and global proper- 
ties. We have been particularly interested in power-law 
scaling of average properties with the chain size. Local 
properties were on the scale of a bead. We specifically 
studied the inter-bead separations, angles, shape and po- 
sitional dynamics. Global or macroscopic properties were 
on the scale of the entire chain. In particular our study 
probed the size and overall motion of vibrated rings. 

The size of various parts of the experimental appara- 
tus limits our study to chains of moderate bead numbers 
typically N < 300 for rings. Consequently the rings are 
analogous to relatively short self-avoiding walks. Despite 
the relatively short chain lengths, we are still able to de- 
duce reasonable scaling laws. The mean global properties 
typically show power-law scaling over the modest range 
in N with cross-over effects at small N and pronounced 
finite-size effects at large N. Generally in between the 
extremes, power-law scaling appears to provide an ade- 
quate description. 

The system of vibrating granular chains presents a 
fertile ground for experimentally studying nonequilib- 
rium statistical mechanics. It is particularly well-suited 
since it permits simultaneous measurements of micro- 
and macro-scopic variables. Our future experiments will 
primarily concern this aspect in the setting of multi-chain 
experiments. We expect to study groups of rings and 
free-end chains and measure the average individual and 
collective properties as they vary with density and con- 
formation. The variation of the average individual or 
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microscopic properties with the macroscopic density wih 
elucidate the statistical mechanics of this system. Rings 
of different A'' have different internal degrees of freedom. 
A tight ring with = 8 is almost circular and has small 
deviations from this configuration under vibrational ex- 
citation. Its space of configurations is small and we may 
think of it as having almost zero internal degrees of free- 
dom. Rings of larger A^, however, have shapes that de- 
part significantly from simple conformations and conse- 
quently a large configuration space and thus many inter- 
nal degrees of freedom. The richness bestowed by the 



ability to vary these internal degrees of freedom will per- 
mit the detailed study of their bearing on macroscopic 
properties of a collection of such rings. Other cooperative 
behavior arising from rings of different sizes, a collection 
free-end chains and self-interactions between different re- 
gions of very long chains afford abundant possibilities for 
studying nonequilibrium statistical physics. 

We thank Mike Rivera and Brent Daniel for useful sug- 
gestions regarding the experiment and Matt Hastings. 
This research is supported by the U.S. DOE (W-7405- 
ENG-3 
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